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Abstract 

Two compressible immiscible fluids in ID and in the isentropic approximation are considered. 

The first fluid is surrounded and in contact with the second one. As the Mach number of the 
first fluid vanishes, we prove the rigorous convergence for the fully non-linear compressible 
to incompressible limit of the coupled dynamics of the two fluids. A key role is played by 
a suitably refined wave front tracking algorithm, which yields precise BV, L 1 and weak* 
convergence estimates, either uniform or explicitly dependent on the Mach number. 
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1 Introduction 

This paper is devoted to the compressible to incompressible limit in the equations of isentropic gas 
dynamics, a widely studied subject in the literature, see for instance the well known results nans 
nuns, the more recent [ TB] , the review m with the references therein and the monograph ma 
for the Navier Stokes equations. For Euler equations, the usual setting considers regular solutions, 
whose existence is proved only for a finite time, to the compressible equations in 2 or 3 space 
dimensions. As the Mach number vanishes, these solutions are proved to converge to the solutions 
to the incompressible Euler equations. 

Consider for instance the isentropic Euler equations in the three dimensional space: 

f dtp + V • {pu) = 0 p( p ) > 0 , P'(p)> 0, 

\d t (pu) + V- (pu®u) + VP(p) =0, (t, x) € [0,+oo[xR 3 . 

where p is the fluid density, u is its speed and P(p) is the pressure. For smooth solutions, this 
system is equivalent to 

(d t p + u • Vp + pV • u = 0 
[d t u + u-Vu+±\7P(p)= 0. 

The Mach number is the ratio between the speed of the particles and the sound speed; it can be 
introduced into the equations in at least two different ways US- 

First, following M, since the incompressible limit can be understood as the limit when the 
Mach number tends to zero, one begins by rescaling the fluid velocity u, —>■ ku where n is a 
small parameter that eventually converges to zero. In order to capture the motion of the particles 
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traveling with a small speed of order of re one needs a space-time rescaling, j —> re f , which allows 
to obtain, in the rescaled variables, the system 


{ dtp + u ■ V p + pV ■ u = 0 
d t u + u ■ Vu + idrVP(p) = 0 . 


( 1 . 2 ) 


Alternatively, the same system is considered in m, but motivated by the following approach, 
see nun]. Consider fluids having equations of state P K (p), parametrized by re, such that the 

speed of sound JP' K {p) —> +oo as re —> 0: 


{ dtp + u ■ Vp + pV • u = 0 
d t u + u-Vu+ ±VP« (p) = 0. 


(1.3) 


The two approaches coincide if the one parameter family of pressure laws P K ( p ) satisfies 

p'Ap) = \p\p), ( 1 - 4 ) 

rZ 

where P is the fixed pressure law as in m- 

In the incompressible limit, the density is constant in time and space so that the functional 
dependence of the pressure on the density is lost. Therefore, it is convenient to use the pressure 
instead of the density as unknown variable. Since P K (p) > 0, we can take the inverse function 

R k (p) = (^Pk'J (p) and rewrite (11.31) using the pressure p as unknown: 

I feify [ dtP + u ■ Vp] + v ■ u = o 

[d t u + u ■ Vu + Tj^yVp = 0 . 

As re —>■ 0, P' K (p) —> +oo, therefore R' K (p) —► 0, and R K (p) —> p, where p is the constant density 
at the incompressible limit. Formally, we get the incompressible equations 

! V -u = 0 

d t u + u ■ Vu + 4 Vp = 0 . 

In [T2] H3] this limit is proved to hold for smooth solutions and small times. The heart of the 
matter is finding energy estimates, uniform in the small parameter re. 

Here, we obtain similar convergence results, in a ID setting, for all times and within the 
framework of merely BV weak entropy solutions. 

The next section describes the physical setting. Section [3] presents the key estimates and the 
main convergence results. All technical details are deferred to Section [IJ 

2 Two Immiscible Fluids 

In a ID setting, an incompressible fluid behaves like a solid and its speed is constant in space. 
Therefore, we consider two compressible immiscible fluids and let one of the two become incom¬ 
pressible, yielding a singular limit for a free boundary problem. Below, we consider a volume of a 
compressible inviscid fluid, say the liquid, that fills the segment [a(t),b(t)} and is surrounded by 
another compressible fluid, say the gas, filling the rest of the real line (see Figure [T]). We assume 
that the gas obeys a fixed pressure law P g (p), while for the liquid we assume a one parameter 
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Figure 1: The two immiscible fluids: the liquid is in the middle, while the gas fills the two sides. 

family of pressure laws P K ( p) such that P' K ( p ) —> +oo as re -» 0. The total mass of the liquid is 
fixed: fa(t) p (t, x) dx = to. Since the two fluids are immiscible, the introduction of the Lagrangian 
coordinate z and of the specific volume r is a natural choice EZI: 


(t,x) = [ d£, 

J a(t ) 


r=\ P g {r) = P g (^j, P K (r) = P B ^ ) . (2.1) 


In these coordinates, the liquid and gas phases become the fixed sets (see Figure O 

£ = ]0,to[ and <5 = M\]0 ,to[. 
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Figure 2: In Lagrangian coordinates, the boundaries separating the two fluids are fixed. 

On P g (r) and P K (t) , we require the usual hypotheses and the incompressible limit assumption: 


P g (r), P K (r) > 0; PL (r), P' K (r) < 0; P'' (r), P" (r) > 0; P' (r) 




->• —oo. 


( 2 . 2 ) 


■ g \ " J i k \ ' j "" yj i g J i k \ ' J ^ ^ g 

In the isentropic approximation, the dynamics of the two fluids is described by the p-system [S] 

jdtr — d s v — 0 where P.( Z , X) = / P ”M f “« £ (2.3) 

\d t v + d z P K (z,t) = 0, \-fff( r ) for z e Q, 

v{t, z ) being the speed of the fluids at time t and at the Lagrangian coordinate 2 . 

The Rankine-Hugoniot conditions for CD , applied at z = 0 and z = to, imply the following 
interface conditions (conservation of mass and momentum) for a.e. t > 0: 

j V (t , 0—) = V (t, 0 + ) j V (t , TO—) = U (t, TO+) 

[Pg (r(t, 0-)) = P K (t (t, 0T)) , \P K (r (t, to—)) = P g (r{t,m+)) . 

In other words, the pressure and the velocity have to be continuous across the interfaces. Hence, 
the pressure is a natural choice as unknown, rather than the specific volume. Therefore, we 
introduce the inverse functions of the pressure laws 


Tgip) = Pg ip ), TM = P* ip ), Xip)^Mo, 
the last limit being a consequence of (12.2D . Rewrite system (12.31) with (p,v) as unknowns 

<a,TA*.P)-a,v = o where r>(2iP) = /r.(ri f»t Z£ £ 

1 + a 2 p = 0 , { 7^ (j>) for z e G . 


(2.4) 


(2.5) 


3 


















The conditions at the interfaces become continuity requirements on the unknown functions: 


r«(t,0-) = „(f,0+) (v ((, m—) = v (t, m+) forae(£0 (26) 

\p (t, 0 —) = p (t, 0 +) [p ( t , iTi ) = p (i, m+) 

The choice of these unknowns significantly simplifies the study of the Riemann problem at the 
interfaces. 


Particular care is necessary to select the one parameter family of pressure laws, the main 
constraint being the validity of ED for all k. Indeed, ED ensures that we recover the same 
equations obtained through scaling and studied in Q21 [13]. The family P K ( p ) = ^P(p) chosen 
in m diverges to +oo as k —> 0. This is not a problem when studying only one fluid as in m 
because the pressure enters the equations only through its gradient. In our case, the value of the 
pressure is very relevant, since it enters the interface conditions (12.61) . Therefore, we cannot allow 
the pressure to grow nonphysically to +oo. We fix the density p of the incompressible fluid in the 
limit and impose that the pressure at that particular density p is a constant, independent of k: 


P K (p)=p , for all k € ]0,1[ . 


(2.7) 


Now, choose a fixed pressure law P = P(p) (for instance, an admissible choice is the usual y-law 
P (p) = k p 1 with 7 > 1 ) and apply conditions ( 11 . 41 ) and ( 12 . 71 ) to get the following expression for 
P K (p), with P{p) =p: 

P K (p) = P + \ \P (p) - P 
K Z L 

which, with the substitution p = -, becomes: 


? 

( 2 . 8 ) 


(2.9) 

= P -1 , we have 


= T(p) = i=T. 

( 2 . 10 ) 


In \ 2 > U), (12.101) is approximated linearly: 

T [p + k 2 (p - p)) w T (p) + k 2 T' (p) (p~p) = f+ k 2 T' (p) (p-p) , 


( 2 . 11 ) 


so that the liquid phase turns out to be governed by a linear system. This approximation makes 
all the estimates simpler. Here we study the Cauchy problem in the fully non linear case 


d t T K ( z,p ) - d z v = 0 
d t v + d z p = 0 , 


where T K {z,p) 


T{p + k 2 (p — p)) for z £ C 
Tg (p) for z eg. 


( 2 . 12 ) 


Colombo and Schleper in [83 Theorem 2.5] proved that for any fixed small n > 0, there exists a 
Lipschitz semigroup of solutions to (12.121 ) . but their estimates are not uniform with respect to n. 
Therefore, as k —> 0 the Lipschitz constant of the semigroup could blow up and its domain could 
shrink, becoming trivial. Here, we provide a full set of new estimates either uniform in k, or with 
the dependence on k made explicit. To this aim, we substantially improve the wave front tracking 
construction in m is], devising and exploiting a different parametrization of the Lax curves. 

The main result of this paper states the rigorous convergence at the incompressible limit in 
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the liquid phase of the solutions to (12.121) to solutions to 


dtTgip) -d z v = 0 
dtv + d z p = 0 




gas; 


< £ _ p(t,0-)-p(t,m+) 

m 


liquid; 


v (t, 0—) = v(t) 

v(t,m+) = v{t) 


immiscibility and mass 
conservation. 


(2.13) 


Note that the liquid speed v(t) is independent of the Lagrangian variable z. In this very singular 
limit, the sound speed in the liquid phase tends to +oo; the density converges to a fixed reference 
value p- the graph of the pressure law P K (r) becomes vertical and the eigenvectors of the Jacobian 
of the flow tend to coalesce. Moreover, the pressure in the liquid wildly oscillates but, remarkably, 
we are able to prove the weak* convergence of the pressure to the linear interpolation of the 
traces of the pressure at the sides of the liquid region, as is to be expected based on physical 
considerations. A linear example, where all the components of this singular limit can be explicitly 
computed, can be found in j5j. 

Recall that problem (12.121) . respectively (12.131) . is well posed in L 1 globally in time, see [HI 
Theorem 2.5], respectively ]3j Theorem 3.6]. 


3 Main Result 

Throughout, we require that the pressure law P g in the gas phase and the one parameter family 
of pressure laws P K in the liquid phase, as defined in (12.91) . all satisfy the condition 

(P): P e C 3 (]0,+oo[;]0,+oo[), P' < 0 and P" > 0. 

The standard choice p{r) = fc/r 7 satisfies this condition for all k > 0 and 7 > 1. 

As a starting point, we provide the rigorous definition of solutions to (12.51) . with reference to |H 
Chapter 4, Definition 4.3 and Admissibility Condition 2]. 

Definition 3.1. Fix T > 0 and k > 0. By weak solution to (12.51) we mean a map 

(p, v) G C° ([0, T]; (L, 1 ^ n BV)(R; R+ x R)) 

such that m holds in distributional sense. The weak solution u is a weak entropy solution 
to (12.51) if both its restrictions to C and to Q are weak entropy solutions in the sense of 
Definition 4-3]. 

Introduce the mathematical entropy flow q = pv of m , the equalities (p, v)(t, 0 —) = (p, v)(t, 0 +) 
and (p, v)(t,m —) = (p,v)(t,m+) (consequences of the Rankine-Hugoniot conditions) imply that 
the entropy flow is continuous and hence that the entropy is conserved across both interfaces. 

In the case of (12.131) . we recall [5J Definition 2.5]. 

Definition 3.2. Fix T > 0. By a solution to (12.131) we mean a pair of maps 

(p*,«*) e c°([ 0 ,r];(L( oc nBV)(g ; i + xi)) 

Vl € W 1,oo ([0, T]; R) 


such that: 


1. (p*,v*) is a weak entropy solution to 


dtTgip) -d z v = 0 
dtv + d z p = 0 


in [0, T] x Q; 
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2. for a.e. t£ [0,T], v t (t) = A (p*{t, 0-) - p*(t, m+)) ; 

3. for a.e. t £ [0, T], v*(t, 0—) = v*(t , m +) = vi(t). 

The existence of solutions to (12.1211 follows from the next theorem, that also provides the basic 
estimates for the subsequent compressible to incompressible limit. In this context, a natural re¬ 
quirement is the smallness of the total variation of the initial datum. Aiming at the incompressible 
limit, it is natural to introduce the weighted total variation 

WTV k (p, v) = TV (p, R) + TV («, Q) + — TV (v, C) (3.1) 

K 

whose boundedness requires that the initial total variation of the particles speed in the liquid 
vanishes with n. 

Theorem 3.3. Fix the total mass of the liquid m > 0 and a pressure p 0 > 0. Let P, P g satisfy (P), 
define T g as in (12. 411 and T K as in (12.101) . Then, there exist positive S , A, L, k* with ac* < 1 such 
that for any n £ ]0,ac*[, for any initial datum (p,v) £ L 1 1 oc (R;M + x R), under the assumptions 

WTV K (p,«) < 6, ||p-po|| L =o( R;R ) < s, v(0+)=v(0) and v(m~) = v{m) , (3.2) 

problem ( 12 . 121 ) with initial datum (p,v) admits a weak entropy solution (p K ,u K ) in the sense of Def¬ 
inition^^ defined for all t £ R + . Moreover, since the specific volume is r K (t,z) = (z,p K (t, z)), 
for any t,t\,t 2 > 0 

WTV k (u K (f, •)) < A, 

TV ( P *(t, ■),£) < A, 

TV (v K {t,-),C) < kA, 

TV (r K (f, •), C) <« 2 A, 

TV(p*(i, •),£?) < A, 

T V(v*(t,-),g) < A, 

T V(^(t,-),g) < A, 

for any z, z\, z 2 £ C 

T V(p*(;z),R+)< A, 

TV (v K (-,z),R+) < A, 

TV (t k (-,z),1+)<kA, 


Sc 

Ic 

Ic 

U 

/a 

lo 


\p K {t 2 ,z) -p li (t 1 ,z)\dz < ip 1 1 2 ~ h \, 
\v K (t 2 ,z) - v Kj (ti,z)\dz < L\t 2 -ti\, 
\r K (t 2 ,z) - r K (ti, z)\ dz < kL \t 2 - ti|, 
\p K (t 2 ,z) -p K {h,z)\dz < L\t 2 -ti\, 
| v K (t 2 ,z) - v K (t 1 ,z)\dz < L \t 2 — h \, 
|r K (t 2 ,2)-r K (ti,x)|d2< Tv |t 2 — ii|; 


(3.3) 


J R + \v K {t,z 2 )-p K {t,z 1 )\ dt < L \z 2 - Zi \, 
f R + \v K ‘{t.,z 2 )-v K {t,z 1 )\dt< kL\z 2 -z 1 \, 
f R + \T K (t,Z 2 ) - T K (t,zi)\dt<K 2 L \z 2 - Zi \; 


(3.4) 


for any z, Z\,z 2 £ Q 

TV (p K (-,z),R+) <A, 

TV (v K (-,z),R + ) < A, 
TV (t k (-,z),R+) < A, 
for any z, z\, z 2 £ R 

TV (p"(-, *),»+)<£, 

TV (v K (-,z), R+) < A, 


f R + \p K (t,z 2 ) -p K (t, 3 i)| dt <L \z 2 - zi\ , 
f K + | v K (t, z 2 ) - v K (t, zi)\dt <L\z 2 - zi|, 
J R + I T K (t, z 2 ) - r K (f, Z\ ) | d t<L \z 2 - Zi\ ; 

f«+ \p K (t,Z2) ~P K (t,z 1 )\ d t<± \z 2 - 2 i| , 
f R + \v K (t,z 2 ) - v K {t,z i)| dt< L\z 2 - zi\ . 


(3.5) 


(3.6) 
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The above existence result can be completed with uniqueness and Lipschitz continuous depen¬ 
dence of the solutions on the data exploiting the results in [ 8 ] Theorem 2.5]. Note however that 
the estimates provided therein, differently from the ones presented here, are not uniform in ft. 

We now pass to the key limit k —> 0. 

Theorem 3.4. Fix the total mass of the liquid m > 0 and a pressure p 0 > 0. Let P, P g satisfy (P), 
define T g as in (12.41) and T K as in (12.101) . Let S, A. L and k * be as in Theorem 1 3. ,‘A For any 
v 0 £ R and (p,v) £ L 1 1 oc (R;R + x R) satisfying 


IIp-Po|| l =o (R;R) < 5, TV (p) + TV (v) < 6 and v(z) = v a Vz€[0,m], (3.7) 


the Cauchy problem for ( 12 . 121 ) with initial datum ( p, v ) admits for any ft £ ] 0 , k*[ a weak entropy 
solution (p K ,v K ) satisfying ( 13 . 31 ) ( 13 . 41 ) ( 13 . 51 ) ( 13 . 61 ) . 

Moreover, there exist functions 

p* £ C°(K+;(L 1 1 oc nBV)(g;l+)), Pl £ L°°(R+ x £;R+), 

v* £ C° (R + ;(L 1 1 oc nBV)(l;l+)) , Vi £ W 1 '°°(R+;R), 

such that (p*,v*g) and vi solve ( 12 . 131 ) with initial datum 

(p*,v*)(0,z) = (p,v)(z) a.e. z £ G 

vi{ 0 ) = v 0 

in the sense of Definition \3.2l Up to subsequences, as k —> 0, 
v K (t,-) —> v* (t, ■) in Iq oc (R;R), t> 0 

v K (-,z) -> v*(-,z) in L 1 1 oc (R+;R), 

TV (v*(t, -),R) < A, f R \v*(t 2 ,z) - v*(h,z)\dz < L\t 2 -h\, t,t 1 ,t 2 >0, (3.8) 

TV z), R + ) < A, f R+ \v*(t, z 2 ) -v*(t, zi)|dt < z,Z!,z 2 £R, 

v* ( t,z) = vi(t), a.e. (t,z) £ R + x C. 


p K (t,-) ->• p*(t,-) in Ll oc (g- R), t> 0 

P K (;z) -»• p*(-,z) in L 1 1 oc (R+;R), z £ Q 

TV -),G) < A, J g \p*(t 2 ,z)-p*(t 1 ,z)\dz < L \t 2 — ti \, t,h,t 2 >0, 

TV (p*(-, z),R + ) < f R + \p*(t,z 2 ) ~P*{t,Z])\dt < L\z 2 -z-l\, z,zi,z 2 £G, (3-9) 

P K ( v) A Pi(;-), in L ao (£xR+;R+), 

Pi{t,z) = (1 - £)p*{t,0-) + ^p*(t,m+), a.e. (t,z) £ R+ x C, 

t k (•, •) —> f, uniformly in C x R + . 

where the specific volume is T K (t,z) = T K [z,p K (t, z)). 


From the Eulerian coordinates’ point of view, the locations of the boundaries of the liquid phase 
can be recovered through a time integration. Let x = a 0 be the initial location of the left interface 
that we keep fixed with respect to n. Since in Theorem 13.41 the initial pressure is chosen inde¬ 
pendently of k, the initial specific volume in the liquid is given by t k (z ) = T [p + k 2 ( p(z ) — p)^j , 
which may depend on k. The total mass m of the liquid is fixed. Hence, the initial location of 
the right interface in general depends on n, say x = b°. Since f K (z) —> t as k —> 0, we have 
bg — > b a = a 0 + mf. Note however that in the particular case of constant initial pressure p(z) = p 
in the liquid, also bf turns out to be independent of k. 

Let a K (f) and b K (t ) be the locations of the interfaces (in Eulerian coordinates) at time t for 
positive «, while aft) and b(t) be the corresponding limits as ft —» 0. Then, we have: 


a K (t) = a 0 + [ v K (5,0) d£ 

a(t) 

= a 0 + [ vi (5) df 


Jo 

ft 


J° 

U 

(3.10) 

b K (t) = b K 0 + u K (5,m) d5 

Jo 

b(t) 

= b a + vi (5) d5 . 

Jo 
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Using Theorem 13.41 we can see that the boundaries of the two phases are Lipschitz continuous 
functions of t. Moreover, as k —> 0, a K —> a and b K —> b uniformly on bounded time intervals. An 
explicit expression for these boundaries and their limit in a linear framework can be found in [5]. 


4 Technical Details 


Throughout, we suppose that P, P g in theorems 13.3113.41 satisfy condition (P) and denote by 0(1) 
a quantity that depends only on P, P g and on uniform bounds on the initial data. 

We define T,T K ,T g as in (12.41) . ( 12 . 101 ) and collect below a few facts about the p-system in 
Lagrangian coordinates using the (p, v) plane. Consider first the gas phase, where 


d t T g {p) -d z v = 0 
d t v + d z p = 0, 


A?(p,v) = 


1 


with eigenvalues 


Vip) 


A 9 2 (p,v) = 

so that the Lax shock and rarefaction curves are, see also [7], 


Vf{p;p 0 ,v 0 ) = 


I 


Vip) 


P<Po 


- ypWdv 

J Po _ 

Vo ~ \j~ ( Tg(p ) - Tg(p o )) ( P - Po) P>Po 

Vo - {Tg(p) -TgiPo)) (p~Po) P<Po 

j Vn V 


P>Po- 


Similarly, in the liquid phase we have 


d t %(p) -d z v = 0 
dtv + d z p = 0 . 


and the Lax curves are 


with eigenvalues 


Viip-,Po,v 0 ) = 


V 2 K (p;p 0 ,v 0 ) = 


v 0 - 


A i(p, v) 

A 2 (P, v) 


K \j T' (p + K 2 (p - p)) 

I / _ 1 

K \ T' (p + K 2 (p ~ p)) 


P <Po 

o - 1 J- {T K (p) - T K {p 0 )) (p - Po) P>Po 

o - \J~ (Tk(p) - T k (Po)) {p - Po) P<Po 

o + / d^ p>p 0 - 


Below we systematically use the parameterizations 

eq ->• Vf{p 0 + op,po , v 0 ) and cr* ->• V*{p 0 + <Ji\p 0 , v a ) 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


of the i- Lax curve, cr, being a pressure difference. Therefore, differently from the usual habit, we 
have that 


Lemma 4.1. Fix L, l with L > l > 0 and n € ]0,1]. The Lax curves (14.41) admit the representation 
Vi(p;Po,v 0 ) = v 0 -K,(p~p 0 )F (U K (p),H K (p 0 )) 

v-iiPiPoiVo) = u 0 + «(p-p 0 )p(n K (p 0 ),n K (p)) 


(4.6) 


8 






















i = 1 

i = 2 

cti < 0 => rarefaction 

02 < 0 => shock 

cti > 0 => shock 

<J 2 > 0 => rarefaction 


Table 1: Types of waves and the signs of the corresponding parameters as in (14. 21) . (14.41) . 



where 


n«(p) 


F(x,y) 


p + k 2 [p - p), 


j \j-V (0X + (1 - 0)y) dd 


V^Vjx) 


x <y, 
x = y, 



—T' (i9x + (1 — $)y) d-d 


x > y. 


(4.7) 


Moreover, 

1. the function F is of class C 1-1 ([Z, L] 2 ; R); 

2. both restrictions F\ x<y an d F\ x>y are °f class C 2 ([Z, X] 2 ; R); 


3. for x,y £ [l,L\, F(x,y) € yJ-V{l) 


The proof follows from standard computations. A property that plays a key role in the sequel is 
that the function F above is independent of n. 

Call F g the function obtained Replacing T with T g in (14.71) . Then, Lemma 14.11 in the case 
k = 1, yields a representation for the Lax curve (ira in the gas phase. 


Riemann Solvers. The wave front tracking algorithm below is, as usual, based on the (possibly, 
approximate) solutions to Riemann problems. 

Throughout, we fix a reference pressure p 0 > 0. By Galileian invariance, in the statements 
below only speed differences will be relevant. 
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Lemma 4.2. There exists a positive 6 such that for all n £ ]0,1] and for any couple of states 


(p,v l ), ( p r ,v r ) with 
( p m ,v m ) satisfying 


P ~Pc 


+ \p r — Po\ < 5 and 


V 1 K (p m ;p l ,v l )=v " 


and 


< kS, there exists a unique state 
V 2 K (p r -,p m ,v m ) = v r . 


Moreover, 


p -p 


\p m -p r \ < 0 ( 1 ) 


P -P 


l V 
v — V 


(4.8) 


A qualitative justification of (14.81) is provided in Figure H) In the liquid region, the Lax curves 
have a slope of order k (see Lemma l4.1[) . hence a jump Av in the velocity generates waves of order 
A v/k. 




Figure 4: Riemann problem in the liquid. Left, in the (t, z ) plane and, right, in the (p, v) plane: 
M + |cr 2 | = 0 ( 1 ) (\p l -p r \ + \v l - v r \/ kJ . 


Proof. Let f = (v r —v l )/n. We apply the Implicit Function Theorem to G(p m ,p l ,p r , £) = 0 where 


G(p m 1 p l 1 p r ,0 = z + (p m - p l ) f (n K (p m ), n K (p z )) - {p r - p m ) f (n K (p m ), n K ( P r )) 

to find p m as a function of (p l ,p r ,f), which is possible since the derivative d p ™G evaluated at 
p m = p 1 =p r = p Q and f = 0 is 

d p mG{p 0 ,p 0 l p 0 , 0) = 2 F (n« (Po) ,n K (p 0 )) 

= 2 \J-T' (p + k 2 {Po - p)) > 2 \J-V (max {p, p Q }) > 0 


Note also that, in a neighborhood of (p 0 ,p 0 ,p 0 , 0), all second derivatives of G are bounded uni¬ 
formly in K, hence the domain of the implicit function contains a neighborhood of (p 0 ,p 0 ,p 0 , 0 ) 

independent of n. Finally, v m can be computed as v m = v l — n ( p m — p l ) F ^II K (p m ), n K (p J )^, by 
Lemma ITU 

Finally, (14.81) follows from G (p l ,p l ,p l , 0^ = G(p r ,p r ,p r , 0) = 0 and the Lipschitz continuity 
of the implicit function. □ 


Note that Lemma 14.21 in the case k = 1 covers the case of Riemann problems in the gas phase 
and slightly improves j4[ Chapter 5]. 

The next Lemma refers to the Riemann problem between the gas, on the left, and the liquid, 
on the right. The symmetric situation is entirely similar. 


Lemma 4.3. There exits a positive S such that for all k £ ]0,1] and for any couple of states 
(p l ,v l ), ( p r ,v r ) with p l — po +\p r — p 0 \ < 6 and v r — v l < 6, there exists a unique state (p m ,v m ) 
satisfying 

Vf{p m -,p l ,v l ) = v m and 


V 2 K (p r ;p m ,v m )=v r . 
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Moreover, 


P -P 


= 0 ( 1) 


P -P 


\p m -p r 1 

= 0 ( 1 ) ( 

P r -P l 

+ 

v l -v r 

-\v m -v r \ 

n 

= 0 ( 1 ) ( 

P r -P l 

+ 

v l -v r 


(4.9) 

(4.10) 

(4.11) 

Av y J 

Proof. Let £ = v r — v l . We apply the Implicit Function Theorem to G(p m ,p l ,p r ,£) = 0 where 

G( p m ,p l ,p r ,0 =f + (p m ~p l )F g (p m ,p l ) -K(p r -p m )F(n K (p m ),n K (p r )) 

to find p m as a function of {p l ,p r ,f), which is possible since the derivative d p ™G evaluated at 
p m = p l = p r = p Q and £ = 0 is 


d p mG{po,p 0 ,p 0 , 0) = \J-Tf(p 0 ) + kF (n« (jpo) ,II K (p 0 )) > > 0. 

Note also that, in a neighborhood of (p 0 ,p 0 ,p 0 , 0), all second derivatives of G are bounded uni¬ 
formly in k, hence the domain of the implicit function contains a neighborhood of (p 0 ,p 0 ,p 0 , 0 ) 
independent of k. Moreover, v m can be computed as v m = v l — (p m —p l ) F g (p m ,p l ), by Lemma l4~Tl 
Concerning the latter estimates, use G(p m ,p l ,p r ,v r — v l ) = 0 to obtain 

m _ , = v l -v r + K (p r — p l )F (H K (p m ),U K (p r )) 

p P Fg(p m ,p l ) + kF (U K/ (p m ),U K (p r )) 


which implies (14.91) and, together with the simple inequality \p m — p r | < p" b — p“ + p“ — p' 

proves (I47TT1 . Finally, the equality v r -v m = k (p r -p m ) F (U K (p m ), U K (p r )), together with { 
proves (14.111) . 


also 


M, 


□ 


Definition of the Algorithm. We modify the standard construction of the wave front tracking 
algorithm, see for instance [J! Chapter 4]. 

First, we identify the state u by means of the pair (p,v). Indeed, we choose to parametrize 
the Lax curves as in (|T 2 P (1441) and, hence, the waves’ sizes are measured through the pressure 
difference a between the two states on the sides of the wave. 

Second, we introduce two strips around the two interfaces z = 0 and z = m, where all 1-waves 
have speed —1 and all 2-waves have speed 1. This, together with [TJ Lemma 2.5], allows to avoid 
the introduction of non-physical waves, significantly simplifying the whole procedure. 

We consider a representative of the initial datum u € (BV DL 1 ) (R,R + x R) such that 
12(0-1-) = 12(0), u(m—) = u(m). Fix e > 0. We approximate the initial datum 12 by a sequence u e 
of piecewise constant initial data with a finite number of discontinuities such that: 

TV(p e ) < TV(p), ||li E -12 || L1 < e, 

TV(l e ;0) < TV(1;0), u e (z) = 12(0) for all z £ \—2e 1 , 2e 2 ], (4.12) 

TV(1 £ ;£) < TV(U;£), u°(z ) = u°(m) for all z G [m — 2e 2 , m + 2s 2 ]. 

Observe that a possible jump at the interfaces 2 = 0 and z = m is assigned to the gas region. 
At each point of jump in the approximate initial datum, we solve the corresponding Riemann 
problem. As usual, see [J Chapter 4], we approximate each rarefaction wave by a rarefaction fan 
consisting of e- wavelets, each with strength less than e and traveling with the characteristic speed 
of the state to its left. On the other hand, each shock wave is assigned its exact Rankine-Hugoniot 
speed. Similarly to what happens in the usual case, there exists a constant S 0 > 0 such that each 
of the above Riemann problems has an approximate solution as long as TV(12) < 5 0 . We introduce 
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two strips around the two interfaces z = 0 and z = m, where all 1 -waves have speed —1 and all 
2 -waves have speed + 1 : 

I~ = [—£ 2 ,£ 2 ] x R + and If = [to — £ 2 , m + £ 2 ] x R + . 

This, together with [fl Lemma 2.5], allows us to avoid the introduction of non-physical waves, 
significantly simplifying the whole procedure. Hence, assign to all 1-waves entering I~ Ul+ speed 
— 1, while all 2-waves entering I~ U If are given speed +1, see Figure0 

Remark that the actual values attained by the approximate solution are not changed, only the 
wave speeds are modified. When exiting these strips, every wave is given back its correct speed. 
By this trick, no interaction among waves of the same family may take place in either of the two 
strips. This construction can be extended up to the first time t\ at which two waves interact, or 
a wave hits one of the interfaces. At time ti, the so constructed approximate solution is piecewise 
constant with a finite number of discontinuities. Any such interaction gives rise to a new Riemann 
problem solved as at time t = 0, if the interaction is in the interior of the two phases, or as 
described in Lemma EM whenever the interaction is along an interface. 

Any rarefaction wave, once arisen, is not further split even if its strength exceeds the threshold £ 
after subsequent interactions, with other waves or with the phase boundaries. The new rarefaction 
waves that may arise at the interfaces are split, if their strength exceeds £, when they exit the strips 
since inside the strips they all travel with the same speed. We can thus iterate the previous 
construction at any subsequent interaction, provided suitable upper bounds on the total variation 
of the approximate solutions are available. As it is usual in this context, see 0 Chapter 7], we 
may assume that no more than 2 waves interact at any interaction point, or that no interaction 
happens at the boundaries of the two strips, thanks to a small modification of the speed of waves 
outside the strips, where necessary. 


Liquid phase 




Figure 5: Left, the strips if and the liquid phase. Right, modification to the usual wave front 
tracking algorithm: waves in the strips I~ and If are assign speed 1, if belonging to the first 
family, and — 1 , if of the second family. 


Interaction Estimates. We recall the classical Glimm interaction estimates, see [4j Chapter 7, 
formulae (7.31)—(7.32)], which hold for any smooth parametrization of the Lax curves: 


<?2 ~ ( a ' + a ") 


erf - (a' + a") 


< 0(1) (JfmT 

< 0 ( 1 ) a'a "| 

< 0 ( 1 ) \a'a"\ 


where we used the notation described in Figure 0 


(Figure 0 left), 

(Figure 0 middle), (4-13) 

(Figure 0 right), 
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Figure 6: Left, an interaction between waves of different families. Center, an interaction between 
waves of the second family. Right, an interaction between waves of the first family. 


Aiming at the convergence result, we need more careful interaction estimates in the liquid 
phase. More precisely, we seek bounds on the constant 0(1) above that allow to control its 
dependence on n. Remark that the choice of parametrizing Lax curves by means of pressure 
differences plays a key role in this improvement. 

Lemma 4.4. There exists a 5 > 0 such that if the interacting waves in Figured hit each other in 
C and all have sizes less than 8, then, the following estimates hold: 


0 Y - a. 


(7 9 (To 


a 2 - ( a ' + <J ") 
r+-(a'+a") + 


< 0 ( 1 ) k 2 

< 0(1) K 2 


< 0 ( 1 ) k 2 \a'c 


a l a 2 


(4.14) 


Proof. Consider first the case of interacting waves of different families, see Figure [HI left. Then, 
with straightforward computations, Lemma 14. 1 1 leads to 

G(a^,a^,af,af) = 0 (4.15) 


where 


Gl(0'i",cr^,cr 1 ,cr 2 ) 
G 2 (cr) t ",cr^,C r f,(xf) 


= + (?2 ~ a l ~ a 2 

= a? F (n K (p 0 + cr) l "),n K (p 0 )j - F (lI K (p 0 + a^),H K (p 0 + erf + <J 2 )) 
-uf F (n K (p 0 + erf + erf ),n K (p D + erf )) + erf F (lI K (p 0 ), U K (p 0 + erf 


Note that by [I] and[2j in Lemma FfTl the function G is of class C 2 and since Ilf (p) = n 2 one can 
compute 


D 2 G(af,aJ,a 1 ,a 2 ) =0(1) 


(4.16) 


Moreover, G(0, 0, 0, 0) = (0,0) and by direct computations, the Jacobian Matrix of G with respect 
to cr+ and a 2 computed at (0,0, 0,0) is 


5 (ff + i(T2+) G(0,0,0,0) = 


yj-r (iUp 0 )) -y/-v (n K ( Po )) 


det 9 (o .+ (7 + ) G(0, 0,0, 0) = 2f/-T' (p + k 2 (p a ~ p)) > 2 \J-V (max{p,p 0 }) > 0 


Hence, the Implicit Function Theorem ensures that (14.151) uniquely defines a map E K of class C 2 
such that (14.151) is equivalent to 

( cr +,cr+) = E K (crf ,crf) 
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for all (a 1 , cr 2 ) i n a neighborhood of (0,0) which can be chosen independently of k. Moreover, 
by (14.161) . 

D 2 E K (<7f,<7^) <C7( 1)k 2 . (4.17) 

By construction, the following equalities are immediate: 

E K (ai,0) = (o-r.O), E K (0,a 2 ) = (0,a 2 ). 

Using (14.171) . compute now 




r 1 


^1 ( CT 1 ! °2 ) a l 

— 

/ 9 Efar^afdtf 
Jo 

°2 


(d a -Z f(trr,^ 2 ") - ^-EftO,^")) 

r / i ^_ CT _E^dv 1 -,^ 2 -)dd / dd 

Jo Jo 12 
2 


d ■& 


a l a 2 


= 0(1) K 


'1 u 2 


We now consider the second estimate in (14.141) . corresponding to the case of interacting waves 
both belonging to the second family. With the notation in Figure [ 6 l middle, we have 

G(<ji,ai,<r',<j")= 0 (4.18) 

where now 

Gi(crf , erf , a', a") = erf + erf - er' - a" , 

G 2 (a^ 1 cr^,a',a") = erf F (ll K (p 0 + af),H K (p 0 )^j -erf F (ll K (p 0 + a^),U K (p 0 + + crf)) 

+a'F (n K (p 0 ), n K (p 0 + a')) + a"F (II K (p a + a'),U K (p 0 + er' + a")) . 

Note that by [D [ 2 ] in Lemma S3 the function G is of class C 2 and since II(. (p) = n 2 one can 
compute again 

D 2 G(at,at,a',a") = 0(1) k 2 . (4.19) 

Moreover, G(0,0,0, 0) = (0, 0) and by direct computations, the Jacobian Matrix of G with respect 
to erf and erf computed at (0,0, 0,0) is, as before, 


D 


(.+ , 4 ) 0 ( 0 , 0 , 0 , 0 ) = 


-v (n M) -J-r (n K (p D )) 


Hence, as before, the Implicit Function Theorem ensures that (14.181) uniquely defines a map E* 
of class C 2 such that (14.181) is equivalent to 

(o-f ,crf) = E K (er',er") 

for all (er', er") in a neighborhood of (0, 0) which can be chosen independently of k. Moreover. 


D 2 Y, K (a',a") < Q{ 1) k 2 . 

L°° 

By construction, the following equalities are immediate: 

E K (<r / , 0) = (0, er'), E K (0, a") = (0, er "). 


(4.20) 
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so that, using (14.201) 


E'»")-(0,a / + <7")|| 

(E>', a") - (0, a' + a")) - (E0) - (0, a')) | 

0 -"E k (<7 / , ’da") — (0,1)) dtf 

d^E K (a', da") - (0,1)) - (3^E K (0, tfa") - (0,1))) dtf 
T, K (d'a',da")dddd' 

= 0(1) k 2 \a' a"\ , 

completing the proof of the second estimate in (14.141) . The case of two interacting waves both 
belonging to the first family in Figure [Gj right, is entirely similar. □ 

The estimates on the waves’ sizes in the case of interactions involving the interfaces are as 
follows. 




Lemma 4.5. There exist positive 6, c and k* < 1 such that, if all the interacting waves in Figure [7| 
have strength less than S, then the following estimates hold: 



< 0(1) K 

°4 

+ (1 + 0(1)(k + 6)) 

CT 2 

7 2 

< (1 — ck) 

°4 

+ (2+ 0(1) 5) 

7 2 


uniformly for all k £ ]0, «»[. Moreover: 


<?1 + ° 2 ~ = + °2 


and 


a 1 = 0 


0 - 2 =° 


r 



-f- 



0 =► 

CT 2 

.+ ' 

°1 


°2 


’2 


’l ^ 


cr 1 a{ > 0 and a 1 aj > 0 


(4.21) 


(4.22) 


Proof. In the present case, we have 


G(a 1 ,^ 2 ^! ,&2 ) = 0 


(4.23) 


where 


Gi(a^,a^,a 1 ,a 2 ) = af + alf - a 1 - a 2 

G 2 (a+,a+,af,af) = af F g (p 0 ,p 0 + af) - Kaf F (ll K (p 0 + af + af),H K (p 0 + af)^j 

+<?t F g (p 0 + af,p 0 ) - na% F (lI K (p 0 + cr^I^^o + af + a£)^ 
with G(0,0, 0, 0) = 0 and Jacobian matrix 


9 KiCT2 +) G(0, 0,0,0) = 


yJ-V(po) -Kyj-V (n K ( Po )) 
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Figure 7: Notation for the proof of Lemma 14.51 o(~, coming from the liquid phase, and o 2 , coming 
from the gas phase, hit against the phase boundary generating er^ in the liquid phase and af in 
the gas phase. 


det 3( ct + jCT +)G(°, 0 , 0 , 0 ) = yJ-T'iPo) + Ky /-V (n K (p 0 )) > tJ-T^Po) > 0 


and moreover 


D~G(a{,aJ,a 1 ,o 2 ) = 0(1) 


uniformly in re, which allows to apply the Implicit Function Theorem in the same neighborhood 
of radius 5 for all small re, yielding a map £ k (°Tj o 2 ) = (cr suc h that 


D 2 E K (o^ 1 o 2 ) = 0( 1) 

locally in (cr^,a 2 ) and uniformly in re. Moreover, 

D£ K (0,0) 

^+,4)^(0, o,o,o)] _1 D^-^Gi 0,0,0,0) 

-v(n K (p 0 )) i 

-i 


(4.24) 


-T'(po) 


kJ-v{um) -J-T'(p 0 ) 


\/~T^{p 0 )+Kyj-T' (n re (po)) 


y/-r^po)+^-r'(nuPo)) 


which shows that the following bound D£ K (0, 0) = 0(1) hold uniformly in re. This, together 
with H4.24p . implies 

0£ K ( t T 1 -, f T 2 -)=0(l) , 


2re^/-T'(n M) -^-T^Po) + KyJ-V ( n «(Po)) 
y / -7;'(p 0 )-re v /-T'(n K (p 0 )) 2yJ-V( Po ) 


so that 


^(a^,o 2 ) = 0(l) 


(4.25) 


since £ K (0,0) = 0. Solve now G^(o^,o 2 ,o x , <r 2 ) = 0 for o+, use the bound (14.251) and the 


estimates 


< S to obtain: 


ov = - 


F g ( Po , Po + cr 2 ) _ F + a i +(J 2 )> n k(Po + )) 


FgiPo + Vl’Po) 


On + 


Fg(p 0 + crt,p 0 ) 


k, a. 


(4.26) 
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(4.27) 


F (ll K (p 0 + a^),H K (p 0 + + o 2 


< (1 + ^( 1 ) s) 

= 0(1) K 


Fg(p 0 +af,p 0 ) 
-0(1) K 

(1 + 0(1)( K + S)) 


K (7 9 


■ 0 ( 1 ) K 


which gives the first estimate in (14.211) . To obtain the second one, use G\ (< 7 +, a 2 , o 1 ,o 2 ) = 0 
and (14.2611 (14.271) : 


/ 


1 + K 


F (ll K (p 0 + a^),U K (p 0 + + cr^)) 


Fg(p 0 + crt’Po) 


1 — K 


F (ll K (po + cf 1 +cr 2 ),II K (p 0 + <j 1 )) 


\ 


F g (p o + (rf,Po) 

which implies the second in (14.211) . since for a suitable c > 0 
(n K (p D + af),U K (p 0 +crf + cr£)j 
F g (p 0 + <ri,Po) 


a i + 1 + 


Fg(p 0 ,Po + a 2 ) 

Fg(p 0 + af,p 0 ) 


(4.28) 


(4.29) 


F I 


> c, 


F (ll K (po + 04 + o 2 ),II K (p 0 + 04 )J 


F g (p 0 + crf,p 0 ) 


Fg(p 0 ,Po + cr 2 ) 

F g (p 0 + af,p 0 ) 


> c, 


< 1 + 0(1) S. 


To prove (14.221) . note that in the case a 1 = 0, (14.281) (14.291) imply that <r 2 and er 2 have the 
same sign. On the other hand, by (14.261) (14(271) 


ct ; = -- 


F(p 0l p 0 + cr 2 ) 

F(Po + Vl,Po) 


0 ( 1 ) 


so that a^ and erf have different signs whenever n is sufficiently small, proving the first equality 
on the right in (14.221) . 

Assume now that erf = 0, so that er^ + a 2 = erf. By (14.281) (14.291) . erf and a 2 have the same 
sign for n small. The second inequality in (14.211) then ensures that a 2 < erf and hence also a+ 
has the same sign of erf and a 2 . □ 


Remark that a wave refracted at the phase boundary remains of the same type, whereas the 
reflected wave changes type when it comes from the liquid and remains of the same type when it 
comes from the gas, see Tableland (14.221) . 

At any fixed positive time t, the approximate solution is a piecewise constant function u e (t) = 
X)q X r r - H t is n °f an interaction time, we denote by cr a the size of the wave supported 

at z a and introduce the potentials 
Vg in = E M 

Oi^Gin 


Vg out = E I CTq 

&£Gout 


Vc= E 


aejs? 


Qq — y \ | Oa O 0 1 Q C — y ( | O a Oft | 

(oi,/3)&Ag (a,f3)£Ac 

T = K in Vg in + Vg out + K c Vc + HgQg + K 2 H c Q c , 


(4.30) 
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where Ki n , Kc, Hg and He are constants independent of n to be precisely defined below. Above, 
we denoted 

Qi n 2 -waves supported in ]—oo, 0 [ and 1 -waves supported in ]m, +oo[. 

Gout 1 -waves supported in ]— oo, 0 [ and 2 -waves supported in ]m, +oo[. 

«5? all waves supported in the liquid phase C. 

Ag pairs of approaching waves supported in the gas phase. 

Ac pairs of approaching waves supported in the liquid phase. 

Here, we define as approaching two waves both supported in the same interval ]—oo,0[, ]0,m[ or 
]to,+ oo[, either of the same family and when one of the two is a shock, or of different families 
with the one of the first family on the right. 

Lemma 4.6. There exist weights Ki n , Kc, Hg and He, all greater than 1, k* £ ]0,1 [ and a 
positive S such that, for all k £ ] 0 , k*[ and piecewise constant initial data u e with the corresponding 
approximate solution u e constructed by the algorithm above satisfying Y(u £ (0+)) < <5, the function 
t —> T (u £ (i)) is non increasing. Moreover, calling o a , op the waves interacting at time t and point 
z, with (J a coming from the left, the following estimates hold: 



G G 

AY 

< 

| & Ot & fi | 


= 0 

AY 

< 

~\o a \ - K \op 


e C 

AY 

< 

-K 2 \ o a Op\ 


= TO 

AY 

< 

-K \o a \ - | Op 


(4.31) 


Proof. Denote by C, with C > 1, a positive constant bounding from above all 0(1) appearing 
in (14.131) . (14.141) and (14.211) . Choose 5 > 0 such that 5 < 1/(2C), and u e such that Y (u £ (0+)) < 5. 

Suppose that at time t there is an interaction and that Y ( u e (t —)) < 5. Consider the dif¬ 
ferent interactions separately. Begin with an interaction in Q, as in Figure [Gj using (14.131) and 
definitions (14.301) : 


A Vg in < C o a op 
A Vg out < C o a op 
AVc = 0 


AQg A C |( T a (7^|h | CTq, op | Y 2 | Oo, op 
A Q c = 0 

AY < (CK m + C- \H g )\o a op\ . 


Consider an interaction in the liquid phase, as in Figure [Gj using (14.141) and definitions (14.301) : 


AVg in = 0 A Qg = 0 

AVg„„ t = 0 A Q c < (C k 2 S — l)\o a op <-\\o a op\ 

AVc < Ck 2 \o a op\ AY < k 2 (C Kc — \ He) o a op\ 

Consider now the case z = 0, the case z = m being entirely analogous. By (14.211) . for k + 5 
sufficiently small so that C(k + 5) < 1, it follows, using definitions (14.301) . that: 


A Vg in 

< 

~\o a I 

A Qg 

< 

2d \o a \ + CkS | op | 

A 

< 

2 cr a | + C n\op\ 

A Qc 

< 

3 6\o a \ + 5 \op\ 

ah £ 

< 

3|cr a | — c n \ o p\ 

AY 

< 

[2 - K in + 3 Kc + (2 Hg + 3 k 2 H c )6] |a Q 
+ K \C — cKc + (C Hg + kHc) 6\ \op\ 


To complete the proof, observe that choosing 

1. Kc so that C — cKc < —2; 

2 . Ki n so that 2 — Ki n + 3Kc < — 2 ; 

3. Hg so that C(1 + K in ) — | Hg < —1; 

4. He so that C Kc — \ He < —1; 

5. 6 so that (CHg + H C )S < 1 and (2 Hg + 3 H C )S < 1. 
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ensures that (14.311) holds. The proof is concluded by induction on the interaction times. 


□ 


Lemma 4.7. With the algorithm defined above, if the piecewise constant initial datum u e is chosen 
so that T(it E (0+)) < 6, with 5 as in Lemma \fiU\ (u e (t) being the approximate solution constructed 
above) then there exists no cluster point of interaction points. 

Proof. By contradiction, call f* the first time at which a cluster point (f*, z*) of interaction points 
appears. 

First, assume that z* ^ 0 and z* 7 ^ m. Call U a neighborhood of (t*,z*) not intersecting the 
interfaces z £ {0, to}. The interactions where there are more than one outgoing waves of the same 
family are those where 

• two waves of the same family hit against each other originating a rarefaction fan of the other 
family of total size bigger than e; and 

• a wave hits an interface, resulting in a new reflected rarefaction larger than e which is 
eventually split as it reaches the boundary of the strip. 

Because of the estimates (14.131) . (14.141) . (14.211) and (14.311) . at any of these interactions AT < — ^e. 
Hence, these interactions may take place only a finite number of times. An application of [TJ 
Lemma 2.5] contradicts the existence of (t*,z*). 

Assume now z* = 0, the case z* = to being entirely equivalent. For a small positive 77 , choose 
a trapezoid A/} contained in I~ of the form 


K = 


(f,z)el : t £ ]t* 


77 , t*[ and 



By construction, finitely many waves cross the lower side of A f v and no wave may enter J\f v along the 
two sides. Inside Af v , any wave can generate another wave at most once, when it hits the interface 
z = z*. Inside J\f v waves propagate with speed either 1 or —1 and at interactions between waves 
with different speeds, no new wave is produced. Hence, the total number of interaction points 
inside M v is finite. This contradicts the existence of a cluster point of interaction points. □ 

To ensure that the value of the functional at t = 0+ is sufficiently small in order that all the 
above interaction estimates hold true, we need some conditions on the total variation of the initial 
data. The standard estimates on the solution of the Riemann problem (see |U Chapter 5]) imply 
that, in the gas, it is sufficient that the initial datum has sufficiently small total variation. On the 
other hand, in the liquid, the estimates on the Riemann problem depend on the small parameter 
k, as shown in dm , see also Figure [4] All this justifies the introduction of the weighted total 
variation m- 

Lemma 4.8. Consider 5 as defined in Lemma \).S\ and let (p e , v e ) = u e : ffi. — > R + xl be piecewise 
constant, continuous at z = 0 and z = to such that || p e — p 0 || L oo < U uE the approximate 
solution constructed above, then, there exists a positive C, which can be chosen bounding from 
above all 0(1) appearing in (14.131) . (14.141) and (14.211) . such that 

i WTV k (to) < Y(u £ (0+)) < C WTV K (u £ ) 

with WTV k defined as in IO and T as in (14.301) with the weights Ki n ,Kc,Hg and He chosen 
as in Lemma Im 

Proof. Let <r Q be the sizes of the waves in u e (0+) and z a be their locations. Consider the estimate 
on the left. The strength of a wave is the absolute value of the pressure difference between the 
states on its sides, therefore, because of the weights’ choice in Lemma 14.61 (they are all greater 
than 1), we have 

TV(p £ ) < TV(p £ (0+)) = Y, M < T(u £ (0+)). 
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The slopes of Lax curves in the gas do not depend on k (Oil . hence, along a Lax curve, the jump 
in the speed is uniformly controlled by the jump in the pressure: 

TV (v E ] g) < TV(u £ (0+); G) < 0(1) ^ |a«| < 0(1) T(u £ (0+)). 

a 

Finally, in the liquid we use (14.61) which shows that along a Lax curve in the liquid, the jump in 
the speed is controlled by k times the jump in the pressure: 

TV(u £ ;£) < TV(u £ (0+);£) < 0(1) K^\a a \ < 0(1) kT(u £ (0+)) . 


This concludes the proof of the left estimate. 

Passing to the right inequality, recall the usual bound T(u £ (0+)) < 0(1) ]E Q |c Q | which clearly 
holds also for T as defined in (14.401) . Proceed using the classical estimate for the solutions to the 
Riemann problems in the gas and (14.81) in the liquid: 


T(u £ (0+)) = 


0(1) J2\ 

ot 

/ 


E m + EI 


= 0 ( 1 ) 

= °( l ) ( ^2 (IP" ( z <x + ) ~ P £ (Za~)\ + \v e (Z a +) - V E (z a -)\^) 


, z a eg 


z a e£ 

= 0(1) WTV K (u E ), 


+ ^2 (j^^T) -p e (z a -)\ + ^|6 e (Za + ) - V S (z a -)\j 


completing the proof. □ 

Proposition 4.9. Fix a positive pressure p 0 and let P 9 ,P satisfy (P). There exist constants 
<5, A, L, k* > 0, with k* < 1, such that, for any k £ ]0,k*[, for any piecewise constant ini¬ 
tial datum u = (p,v), continuous at the points z = 0, z = m, satisfying WTV k (m) < <5 and 
||p —p 0 || L oo < 6, the wave front tracking approximate solution zt K,£ = (p K ’ s ,v K ’ e ) to the Cauchy 
problem for (12.121) can be constructed for all times t > 0. Moreover, given the specific volume as 
T K,e (t, z) = T k (z,p K,e (t, z)), the following estimates hold. 

For any t, ti,t 2 > 0 


wtv k 

(u^(t,-)) 

< 

A, 












TV ( p K - 

e (t,-),c) 

< 

A, 

Sc 

\p K 

' e (t2, 

z) - 

-p K - 

' E (h, 

z)\ 

dz 

< 


- ^1, 

TV ( v K > 

E (t,-),c) 

< 

kA , 

Ic 

\v K 

’ e (h, 

,z) - 

- v K 


■z)\ 

dz 

< 

L \t2 - 

- ill, 

TV (t k - 

e (t,-)X) 

< 

k 2 A, 

Ic 

\t k 

,e (ti- 

>z) ~ 

-t k 

,£ (ti, 

,z)\ 

dz 

< 

kL |t 2 - 

- ii|, 

TV (p K ’ 

e (tr),G) 

< 

A, 

U 

\p K ■■ 

E (h, 

*)- 

-P K ’ 

e (ti, 

z)\ 

dz 

< 

L |t:2 - 

- ii|, 

TV ( v K • 

e (t,-),G) 

< 

A, 

S 0 

\v K - 

,e (t2, 

*)- 

- V K - 

• e (ti, 

z)\ 

dz 

< 

L |t 2 - 

- ill , 

TV (r K - 

e (t,-),G) 

< 

A, 

Is 

| t k 

' E (h, 

z) - 

- t k 

,E (ti- 

a) 

dz 

< 

L 1f 2 - 

-ill • 
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(4.32) 







For any z £ £, z\, Z2 € C \ ([—e 2 , e 2 ] U 

[m — e 2 ,m + e 2 ]) 

TV (p*’ E (-, 

z),R + )< f, 

f R + jp K,£ (t,z 2 ) -p K ’ £ (t,z 1 )j d t < L\z 2 - zi 

TV (v K ’ e (-, 

<i 

VI 

m 

w 

f R+ \v K ' e (t, z 2 ) - v K ’ e (t,zi)\ d t < kL\z 2 - zi 

TV (r K - E (-. 

, z), K. + ) < kA, 

f R + \r K ’ £ (t, z 2 ) - r K ’ £ (t,zi)\ dt < k 2 L \z 2 - zi 

For any z, z\, z 2 £ Q 



TV (p K ’ E 

(■,«),»+)< A, 

f R + P K,£ (C z 2 ) - p K ’ £ (t, z\)\dt < L\z 2 — zi\ , 

TV (iO- E 

(■,*),R+)<A, 

f R + \v K ’ £ (t, z 2 ) - v K ’ £ (t, zx)\dt <L\z 2 — zi \, 

TV (r K ’ E 

(;z),R+)<A 1 

f R+ \r K ’ e (t , z 2 ) - r K, ’ £ (t 1 zi)\ dt<L \z 2 — zi\. 

For any z,zi,Z2 6 t 



TV (p K ’ E 

{;Z),R+)<$, 

f R+ \p K ' e (: t , Z 2 ) - P*’ £ (t,zi)\dt<±\z 2 - zi \, 

TV (v K ’ E 

(; Z ),R+)< A, 

f R+ I v K ’ £ (t,z 2 ) - v K ’ £ (t,z i)| dt< L \z 2 - Zl \. 


(4.33) 


(4.34) 


(4.35) 


Moreover, the maximal size of rarefaction waves is uniformly bounded by a constant, independent 
of k, times e. 

Proof. Choose 6 as in Lemma 14.61 and k* as in Lemma 14.51 Define 6 = S/C, with C as in 
Lemma 14.81 Using the piecewise constant initial data u, we use the previously described algorithm 
and call w K,E the piecewise constant approximate solution so obtained. By Lemma 14.81 we have 
T (m k,e (0+)) < C ■ WTV k (u) < 5. By Lemma ITTil the map t —> T (u K,E (t)) is not increasing so 
that T (u K ’ e (t)) < 6 for all positive times. Lemma [4771 ensures that u K,£ can be constructed for all 
times t > 0. Again, Lemma 14.81 implies the estimate 


WTV« ( u K ’ E (t .)) < CT (u K ’ e (t,)) < C T (it K,E (0+)) < C 2 WTV k («) < C 2 S = CS. 

The estimates on the total variation of p K,E and v K,£ in (14.321) immediately follow. To obtain 
the bounds on the total variation of the specific volume in the liquid, use (12.101) . The Lipschitz 
continuity estimates in (14.321) are now a standard consequence, see e.g. [4j Section 7.4], since the 
wave propagation speed in the gas is uniformly bounded independently of n and in the liquid (also 
in if) is bounded by 0(1) / k. 

Pass to (14.331) . Observe first that from the proof of Lemma 14.61 it follows that 


6 C 3 5 < 1 and C >1. 


(4.36) 


As usual, we call a a the size of the wave supported at z a . For z £ C, define 


T ,(t) 
w z (t) 

V z (t) 

h(t) 
T (t) 


30 

w z (t) + V z (t) + — T(t) 

K 

Y 

rG[0,£] 

Y m 

a€lz(t) 


or. z a € £ and the wave at z a is of the 


first family and z a > z 
second family and z a < z 


T(«*1t)) 


Note that the sum defining W z is actually a finite sum, since the total number of waves is finite 
by Lemma ITT! We claim that t —> T z (t) is non increasing. Indeed, T z (t) may change its value at 
a time t when: 
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1. A wave with size cr s crosses z and no other interaction occurs. Then, A W z (t) = \Ap(t,z)\ = 
|<r a |, A V z (t) = — |cr a | and AT(t) = 0. Hence, AT z (t) = 0. 

2. An interaction in Q occurs and no wave crosses 0 . Then, AW z (t) = 0, AV z (t) = 0 and 
AT (t) < 0. Hence, AT z (t) < 0. 

3. An interaction in C occurs and no wave crosses z. Then, AW z (t) = 0; calling o a , op the sizes 
of the interacting waves, A V z (t) < Cn 2 \o a cr^ | by (|4.14[) and AT (t) < —K 2 \o a op \ by (14.311) . 
Hence, AT z (t) < C (l — |) K 2 \o a op\ < 0. 

4. A 2-wave with size cr a , coming from Q, and a 1-wave with size op, coming from C , interact at 
z = 0. Then, AW z (t) = 0; bv Lemma l4.51 AV z (t) < (1 — cK)\op\ + (2+C5)\o a \] by Lemma l4.6l 

AT (t) < —|cr a | — «|crja|- Hence, AT z (t) < (1 — ere — 3 C) \op\ + (2 + CS — |cr a | < 0. 

5. Two waves interact at z = m: the same procedure as above applies. 

The remaining times where T z may change value consist in the superposition of two or more of the 
cases considered above and can be dealt superimposing the corresponding inequalities. Therefore, 

TV(p K ’ e (-, 2 )) = sup TV (p K ’ E (-, z); [0, T]) = supWd(T) < sup T Z (T) 

T >0 T>0 T>0 

< T,(0+) = K(0+) + ^T(0+) < f 

provided A > 2 6, completing the proof of the first estimate on the total variation in (14.331) . The 
remaining total variation bounds in (14.331) follow from the estimates 

\Av(t,z)\ < 0(1) k \Ap(t,z)\ and \Ar(t, z)\ < 0(1) k 2 \Ap(t, z)\ 

which hold along Lax curves by Lemma FTTl and (12.101) . The Lipschitz continuity estimates in (14.331) 
are now a standard consequence, see e.g. [H Section 7.4], since the wave propagation speed in 
C \ ([—e 2 , e 2 ] U [m — e 2 , to + £ 2 ]) is of order 1/re. 

The proof of the estimates (14.341) is obtained from that of (14.331) completed above, formally 
setting re = 1 and with obvious modifications to the definition of T z . 

The estimates on all the real line (14.351) are obtained choosing a common upper bound on the 
total variation and a common lower bound on the wave speeds in the liquid, in the gas and in the 
two strips If and observing that z 1 —► u K ’ e ( t , z) is continuous at z = 0, 2 = m for every t > 0 in 
which no wave interacts with the interfaces. Observe that a similar Lipschitz estimate does not 
hold for the specific volume r, since at z = 0 and z = m it is not continuous. 

Finally, the estimate on the maximal size of rarefaction waves follows the lines in [5] Section 7.3, 
Step 5]. Indeed, call o(t) the size at time t of a rarefaction wave in the wave front tracking 
approximation. We claim that, if in the interval [f 0 ,r] the wave does not leave the phase in which 
it is found at time t 0 and does not disappear due to possible interactions with shocks of the same 
family, then |<t(t)| < 6|<r(f 0 )|. 

Indeed, consider the liquid phase, let z(t) be the location of the wave at time t and define 
s (t) = \o(t)\ [l + 6C 2 re 2 H s (f)-t-24C 3 reT(f)] 

= Sae/ s (t) M 

J s (f) = {a: z a £ C, the wave at z a is approaching the wave at zj. 

The function t —> s(t) is non increasing in the interval [t 0 ,T\. Indeed, s(t) may change its value at 
the following times: 
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1. At time t a wave a a interacts with the wave at z(t) and no other interaction occurs. Then, 
by (14.1411 A |(t(t)| < Ck 2 \a(t—)cr a \; A V s (t) = — |<r a |; AT (t) < 0. Hence, by (14.361) . 

As(t) = A |(t(t)| \l + 6 C 2 K 2 V s (t+) + 24C 3 nT(t+)] + |d(f-)| \&C 2 K 2 AV s {t) + 24C 3 KAT(t) 


< Ck 2 \<r(t— )cr Q 

< Ck 2 \a(t— )a a 


1 + 6C 2 k 2 5 + 24C 3 kS — 6C 2 k 2 \cr{t —)| |er Q 
1 + 6 C 2 k 2 S + 24 C 3 k6 - 6 C 


< Ck 2 \cr(t— )a a \ [1 + 1 + 4 — 6 C] < 0 . 


2. At time t, an interaction in Q occurs and no wave crosses z(t). Then, A |cr(t) = 0, A V s (t) = 0 
and AT (t) < 0. Hence, A s(t) < 0. 

3. At time t, an interaction in £ occurs and no wave crosses z(t). Then, A |cr (t) | = 0; calling cr a , 
<jj 3 the sizes of the interacting waves, AV s (t) < Ck 2 | cr a ap | by (14.1411 and AT(f) < —K 2 \(j a ap\ 
by (14.3111 . Hence, 


A s{t) < |cr(t)I 6C 3 k 4 I cr a (T^ I — 24C 3 k 3 Io-qO-^I < |cr(t)| Icrccr^l 6C 3 k 3 (k — 4) < 0 . 


4. At time i, an interaction occurs at z = 0 and no wave crosses z{t). Call a a the size of the 
wave coming from Q, and crp the size of the wave coming from £. Then, A |d(t)| =0. By 
Lemma 14.51 A V s (t) < (1 — cK)\ap\ + (2 + C5)\a a \\ by Lemma PTT7TI AT (t) < — |cr a | — /cjcr^l. 
Hence, 


A s(t) < \<y(t )| 6 C 2 k 2 ^(1 — c/c)|cr^| + (2 + C'<5)|cr a |^ — 24 C 3 k ^|cr a | + K | cr /3|) 

< 6 C 2 k |(r(t)| K ^(1 — CK)\a 0 \ + (2 + C'5)|cr a |^ — 4C ^|cr a | + «|cr/31^ 
<6C 2 n\a(t)\ k | rr^ 3 1 (1 — ck — 4C) + |cr Q | (2k + CkE — 4C) <0 


5. Two waves interact at z = m: the same procedure as above applies. 

The remaining times where s(t) may change value consist in the superposition of two or more of 
the cases considered above and can be dealt superimposing the corresponding inequalities proved 
above. Therefore s(r) < s(t 0 ) which implies 


c ( t )| < W(to) 


1 + 6 C 2 K 2 V s (t 0 ) + 24C 3 K T{t 0 ) 
1 + 6 C 2 k 2 V s {t) + 24C' 3 «T(t) 


< 


1 + &C 2 5 + 24C ,3 (5 


|cr(t 0 )I < 6 |cr(t 0 )| 


This proves the claim in the liquid. In the case of a wave in the gas, the argument is similar: it is 
sufficient to set n = 1 in the definition of s(t) and make the obvious modifications to the map V s . 

Finally, we observe now that when a wave crosses the interfaces, the refracted wave has a 
strength given by the strength of the incoming wave times a constant bounded uniformly with 
respect to k, for instance we can choose 3C (see Lemma [4.511 . Moreover, when a rarefaction is 
born, its strength is less than e and it can cross at most an interface once. Therefore, also the last 
claim of the Proposition is proved with the constant 1944 C 2 e. □ 


Proof of Theorem 13.31 Use 6 , A,L,k* > 0 as defined in Proposition 14.91 and choose any 
k £ ]0, k*[. Fix a suitable sequence £„ strictly decreasing to 0. Approximate the initial datum 
u = (p,v) with an approximate, piecewise constant initial datum u El ’ satisfying (14.121) . so that 
WTV k (u e ") < WTV k (u) < (5, II fF" - Po IIloo < <5. 

Proposition 14.91 ensures that it is possible to construct a wave front tracking £„-approximate 
solution (p K > e ^,t) K > e ^) that satisfies all properties stated therein. 
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Using (14.321) and (14.351) . a repeated application of Hclly Theorem [4J Theorem 2.4], ensures the 
convergence of a suitable subsequence, which we still denote by u Kj,e ‘ / , to a function u K = (p K , v K ) 
in the following sense 


lim II (p K ’ e ",i; 

z^—)-+oo 11 

lim ||(p K ’ e,/ , 

v —»+oo 11 


(p ’ V )(*’ ')|Il 1 ([-M,M];R+xR) 

V ’ z )~(p > v )(■) 2; )|lL 1 ([0,Ar];R+xR) 

(p K ,u K )(0,-) 


0, for any t > 0, M > 0 
0, for any z£l, M > 0 

(P K , ^HO- 


Passing to the limit in (14.321) . (14.331) . (14.341) . (14.351) . we obtain (13.31) . (13.411 . (13.51) and (13.61) . 

Since the bounds on the total variation are uniform in e and since the strength of rarefactions 
is uniformly bounded by a constant times e, standard techniques in wave front tracking pH Section 
7.4] can be used to show that the limit u K is a weak entropy solution to (12.121) in the open regions 
z < 0, 0 < z < m, z > m. By (13.61) . we have that the map 2 —► u K (-,z) is continuous in L 1 , 
in particular it is continuous across z = 0 and z = m. Therefore, u K trivially satisfies there the 
Rankine-Hugoniot conditions and the entropy (in)equality. Hence, u K is a weak entropy solution 
to in all 1+xl. " □ 

Proof of Theorem 13.41 By (13.71) . WTV K (u) <6 so that Theorem 13.31 applies, ensuring the 
existence of a solution u K = (p K ,v K ) to (12.121) satisfying (13.31) . (13.41) . (13.51) and (13.61) . 

Since k < 1, from (IH and m we have for 

TV (v K (t,-),R) < A, f R \v K (t 2 ,z) - v K (t 1 ,z)\dz < L\t 2 -h\, t,t 1 ,t 2 >0, 

TV(u k (-, 2 ;),M + ) < A, J R+ \v K (t,z 2 ) - v K (t,z 1 )\dt < L\z 2 -Z!\, z,Zi,z 2 gR. n 

Helly Theorem ^] Theorem 2.4] implies the existence of a subsequence (that we call again v K ) 
converging to a limit v* in the sense of (13.81) . From the bound in (13.31) on the total variation of 
v K or from the Lipschitz estimate in (13.41) for v K in the liquid, it is straightforward to obtain that 
v*(t,z) = vi(t) for all z £ C and t > 0, where vi(t) is a function which depends on time only, 
completing the proof of (13.81) and of 3. in Definition 13.21 

The same procedure can be carried out for the pressure in the gas region, proving the first four 
lines in 02). Observe that for the pressure, we cannot apply Helly Theorem in the liquid since 
there the estimates blow up as k —> 0. Because of the strong convergence in the gas region of both 
the velocity and the pressure, the limit u* = (p*,v*) satisfies 1. in Definition 13.21 and the initial 
condition u*( 0, z) = u(z) a.e. z € Q. 

The uniform convergence of t k in the liquid is a straightforward consequence of (12.1 HI) and of 
the uniform bound on the L°° norm of p K . 

Since the pressure is uniformly bounded, we have a weak* convergence (possibly passing to 
further subsequences) p K ^p* in L°° (R + xl,E) ]S] Section 4.3 Point C.]. If we define pi = P* R + x£ 
we get the fifth line in (13.91) . 

By (13.31) and (13.61) . the second equation in (12.121) can be written in integral form in [t\,t 2 ] x C: 

pm pm pt2 pt2 

/ v K (ti,z ) dz— v K (t 2 ,z ) d z+ p K (t,0) dt— p K (t,m) di = 0. (4.38) 

J 0 t/0 J t,\ 

Now, we use the strong convergence of both p K and v K in the gas region and the fact that in £, 
v* is constant to obtain 

m [vi (f 2 ) - vi (Zi)] = f p* (t, 0) dt- f p* (t,m) dt. 

Jtx Jtx 

Setting ti = 0 and £2 = t in the last expression above, 

Vi(t)=vi( 0) + — / [p* (s, 0) - p* (s, to)] ds = vi (0) + — [ [p* (s,0-) - p* (s,m+)] ds 
w J 0 171 J 0 


24 
























which means that vi is Lipschitz continuous and satisfies 2. in Definition 13.21 
Observe that the non linear term T K ( z,p K ) converges strongly to 


T *n z \ _ J T for z e C, 

[T g (p*(t,z)) for z £ Q, 

hence we can pass to the limit in (12.121) in distributional sense to obtain 


d t r* - d z v* = 0 
d t v* + d z p* = 0, 


x R. 


(4.39) 


Since in the liquid region v* ( t, z) = vi ( t ) with vi Lipschitz continuous, the second equation in (14.391) 
becomes 

d z p*{t, z) = —vi(t) in R + x C. 

Therefore there exists a measurable function j3(t) such that the function 


Pi{t,z) = —zinit) + P{t) 


can be chosen as a representative of the limit pressure p* restricted to the liquid. This implies the 
existence of the two limits 


lim pi(t, z) = P{t), lirn _pi(t,z) = /3(t) - zinit). 

z —^0 ' z—±m 

The fourth line in (13.91) ensures the existence of the corresponding limits from the gas region: 
lim p*(t, z) = p*{t, 0), lim p*(t, z) = p* ( t, in), a.e. t £ R + , 

2—>-0“ 2—)-m+ 


hence Rankine-Hugoniot conditions for (14.391) applied along z = 0 and z = m imply that the right 
and the left limit of the pressure must coincide along z = 0 and z = m for a.e. t > 0. Therefore, 
we have 


P*{t, 0) = Pit) 
p*it,m) = —mil (t) + Pit) 


for a.e. t > 0, 


which implies the remaining equality to be proved in (13.91) . □ 
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